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We excite ion acoustic waves in ultracold neutral plasmas by imprinting density modulations during plasma 
creation. Laser-induced fluorescence is used to observe the density and velocity perturbations created by the 
waves. The effect of expansion of the plasma on the evolution of the wave amplitude is described by treating 
the wave action as an adiabatic invariant. After accounting for this effect, we determine that the waves are 
weakly damped, but the damping is significantly faster than expected for Landau damping. 



I. INTRODUCTION 

Collective modes in plasmas 1 , are important dynamical 
excitations and can often be used to obtain information 
on plasma density, pressure and temperature. Here, we 
describe the excitation and study of ion acoustic waves 
(IAWs) in ultracold neutral plasmas (UNPs) 2 ~ 4 . 

UNPs are created by photoionizing laser-cooled atoms 
at the ionization threshold. They have ion and elec- 
tron temperatures that are orders of magnitude colder 
than traditional neutral plasmas, so they represent a 
new regime in which to study collective effects, where 
the ions display correlated particle dynamics reflecting 
strong coupling 5-7 . Furthermore, the density distribu- 
tion and temperatures can be precisely controlled and 
probed to enable a broad range of experiments. 

Experiments and theory on collective modes in UNPs 
have explored Langmuir oscillations 8 - excited by RF elec- 
tric fields^. Newly applied techniques of rf absorption 10 
have confirmed the relationship of the Langmuir oscil- 
lations to edge modes and their dependence on non- 
neutrality induced by plasma dynamic s 11 ! 12 . Similar ex- 
periments identified Tonks-Dattner modes in a series of 
RF resonances observed at frequencies above the Lang- 
muir oscillation 13 . A high-frequency electron drift insta- 
bility was observed in an UNP in the presence of crossed 
electric and magnetic fields 1 ^. Numerical simulations 
have identified the possibility of propagating spherically 
symmetric ion density waves 1 ^.. A theoretical treatment 
of IAWs in UNPs, including the effects of strong cou- 
pling, was published^ and ion-acoustic shock waves were 
recently discussed 1 ^. The first experimental excitation of 
IAWs in UNPs was recently reported^. This represents 
the first study of ion wave motion in UNPs. 

To excite IAWs, we create density perturbations by 
spatially modulating the intensity of the laser that pho- 
toionizes the atoms to create the plasma 2 . Doppler- 
sensitive laser-induced fluorescenc e 18 ' 19 is used to study 
the dispersion relation and damping of the waves and to 
measure the velocity distribution of the ions. These tech- 
niques open new areas of plasma dynamics for experimen- 
tal study, including the effects of strong coupling on dis- 
persion relations 2 ^— and non-linear phenomena ^ 15 ' 24 ' 25 
in the ultracold regime. 



Low frequency ion density waves in the absence of a 
magnetic field are described by the familiar dispersion 
relation for frequency ui and wavevector 

/u\ 2 _ k B T e /mj 

\k) ~ l + k 2 X 2 D ( ' 

where rrii is the ion mass, T e is electron temperature, and 
Ad = v £ofcs7e/n e e 2 is the Debye screening length, for 
electron density n e and charge e. We have neglected an 
ion pressure term because the ion temperature satisfies 
Ti <C T e in UNPs. In the long wavelength limit, which is 
the focu s here, th is mode takes the form of an IAW with 
lu = k^/huTe/mi, in which ions provide the inertia and 
electron Debye screening moderates the ion-ion Coulomb 
repulsion that provides the restoring pressure. 

IAWs are highly Landau damped unless Ti -c Te 1 , 
however, they are often observable in high-temperature 
laboratory plasmas 2 ^ - — . Acoustic waves of highly 
charged dust particles, which show similar charac- 
teristics, have been studied experimentall y 27 ' 28 and 
theoretically^ ."^. Beyond fundamental interest, IAWs 
are invoked to explain wave characteristics observed in 
Earth's ionosphere^ 9 - and transport in the solar wind, 
corona, and chromosphere^ .. 

II. EXPERIMENTAL DETAILS 

A. Creation and initial dynamics of an ultracold neutral 
plasma 

Ultracold neutral plasmas are created by photoion- 
ization of laser-cooled strontium atoms in a magneto- 
optical trap (MOT^i - ^ 3 -. Operating on the dipole- 
allowed 1 So — 1 Pi transition of 88 Sr at 461 nm, we rou- 
tinely trap 2 x 10 8 atoms at a temperature of ~ 10 
mK in a spherically symmetric Gaussian density distribu- 
tion, ti (r) = n atoms exp(— r 2 /2cr 2 ), with a ~ 0.6 mm and 
n-atoms ~ 6 x 10 10 cm~ 3 . Prior to photoionization, the 
MOT lasers and magnetic field are turned off, and the 
atoms are allowed to expand to obtain larger samples, 
with corresponding lower densities. 

Photoionization of the atoms is a two-photon process 
performed by two tempo-spatially overlapping, ~10 ns 
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FIG. 1. Overview of UNP dynamics. Note the three distinct 
time scales of electron equilibrium, ion equilibrium and global 
plasma dynamics, (adapted from Ref. 0). 



laser pulses: the first from a pulse amplified cw beam 
tuned to the cooling transition of 461 nm and the second 
one from a pulsed dye laser tuned just above the ioniza- 
tion continuum at ~ 412 nm. The pump for both beams 
is a 10 ns Nd:YAG laser, operating at 355 nm with a 10 
Hz repetition cycle. This process ionizes ^30-70% of the 
atoms, creating an ultracold neutral plasma. For a spa- 
tially uniform photoionization laser intensity, the plasma 
inherits its density distribution from the MOT, resulting 
in initial electron and ion densities (rtn e ~ firjj = no) as 
high as « 4.2 x 10 10 cm~ 3 . Effects of un-ionized atoms on 
the plasma are not significant considering the fast time- 
scales of the experiment and small neutral-ion collision 
cross-sections. 

Due to their relatively small mass the electrons ac- 
quire most of the excess energy from the photoionizing 
beam, while the ions' kinetic energy remains similar to 
that of the neutral atoms in the MOT—. By tuning dye 
laser wavelength, we create electrons with initial kinetic 
between 1 and 1000 K, with a resolution of about IK. 

On a timescale of the inverse electron plasma oscilla- 
tion frequency (l/ui pe ), electrons equilibrate to produce 
a nearly thermal distribution 4 (Fig. [I}, yielding average 
Debye lengths, Ad from 3-30 /jm. Three-body recombina- 
tion can also produce large numbers of Rydberg atoms 
during this initial stage if the electron temperature is 
low enough 3 ^. (We avoid this recombination-dominated 
regime for our study of IAWs.) On a timescale of the in- 
verse ion plasma oscillation frequency disorder- 
induced heatin g 35 ' 36 raises the temperature of the cold 
ions to approximately 1 K. On a longer hydrodynamic 
timescale, given by the ratio of the initial characteristic 
plasma size to thermal velocities, 
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the plasma expands into the surrounding vacuum—. For 
a strontium plasma with a typical initial size of <x(0) = 1 
mm, an initial electron temperature, ^(O) = 50 K, and 
an initial ion temperature, Ti(0) = 1 K, r exp — 14 [is. 

The expansion of the plasma has a large effect on 
IAWs that must be accounted for in any quantitative 
model. Expansion of an unperturbed UNP was stud- 
ied theoretically i n 15 i 35 i 38 ~ — and was seen experimen- 
tally in^ i 32 ' 34 ! 37 ' 41 " — . Fundamentally the expansion is 



driven by the pressure of the electrons coupled to the 
ions through a space charge field. 

Under the assumptions of quasi-neutrality, adiabatic- 
ity, and spherical Gaussian symmetry, the evolution of 
the size of the plasma cloud and the hydrodynamic ex- 
pansion velocity (u exp ) can be written as 



t z 



^(*) = ^(0)(l + -5-) 



7(*) = 



t/r, 



exp 



exp 



u exp {r,t) = 7(i)r, 
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where r is the radial distance from cloud center. The ion 
and electron temperatures (Tj and T e ) drop during the 
expansion due to adiabatic cooling as thermal energy is 
converted into the kinetic energy of the expansion, 
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This expansion is similar to dynamics seen in plasmas 
produced with solid targets, foils, and cluster o 37 ' 44 " — . 



B. Optical diagnostics of density and velocity distributions 

The optical diagnostics used to study UNPs have been 
described in detail previousl y 19 ' . The primary probe 
is Doppler-sensitive laser induced fluorescence (LIE), ex- 
cited by linearly polarized light of frequency v near reso- 
nance, vq, with the primary 2 Si / 2 — 2 Pi/2 transition of the 



Sr ions at A =422 nm. The natural linewidth of this tran- 
sition is 7o/27r = 20 MHz. Wc define x as coaxial with 
the LIF beam, and fluorescence is imaged along the z axis 
onto an intensified ccd camera (Figure [2]) . To minimize 
density variations along the line of sight, the excitation 
light is aligned to the center of the plasma and formed 
into a sheet in the x — y plane with 1/e 2 intensity radius 
w z = 0.625 mm. We derive spatio-temporally resolved 
density and velocity distributions for the ions from the re- 
sulting images, F(x,y,v). Frequency dependence arises 
from the natural linewidth, 70, and Dopplcr-broadening 
of the transition, and images can be taken with time- 
resolution of 50 ns. 

A single fluorescence image can be related to underly- 
ing physical parameters through 



F{u,x,y) = C^ f dz ni (r)^l 
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where the multiplicative factor, C, depends upon col- 
lection solid angle, dipole radiation pattern orientation, 
and detector efficiency; it is calibrated using absorption 
imagin g 19 ' 32 ' 37 . I(r) is the intensity profile of the flu- 
orescence excitation beam, and I sa t is the saturation 
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FIG. 2. Partial experimental schematic showing transmission 
mask along the path of the ionizing beam. The ionizing beam 
is allowed to first pass through the plasma and the transmis- 
sion mask is placed as close as possible to a retro-reflecting 
mirror. Due to optical access constraints the k-vector of the 
I AW is at a 16° angle with the k-vector of the fluorescence 
beam. The fluorescence beam is a sheet perpendicular to the 
imaging axis (adapted from Ref . [49h . 



intensity for linearly polarized light. Taking Clebsch- 
Gordon coefficients for the transition into account, I sa t — 
114mW/cm 2 . Jeff is the sum of the natural linewidth 
and instrumental linewidth. The detuning factor in Eq. 
[7] reflects the Doppler shift (v x /X) of the laser frequency 
for an ion with velocity along the laser of v x . The veloc- 
ity distribution function for the ions, assuming uniform 
ion temperature, is 
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with velocity width s v — yfkBTi/rrii. We allow for a gen- 
eral hydrodynamic velocity u that can vary with position 
and will contain expansion velocity (Eq. [5]) and any other 
ion motion, such as motion induced by an IAW. 

A measure of the ion density is obtained by summing 
the fluorescence signal for a series of images taken at 
equally spaced frequencies covering the entire ion reso- 
nance: 



J dvF{v,x,y) = C^ Jdzm(r)j 



I(r) 



(9) 



For the imaging-laser sheet geometry used here, lu z <C <j, 
one obtains a measurement of the ion density in the plane 
of the laser, 

. . f dv F(u, x, y) , „. 

n(x,y,z = Q)* J , ' (10) 
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Figure [3] shows examples of two-dimensional false color 
plots of the ion density distribution. 



FIG. 3. 2D density profiles of density perturbations in UNPs 
for T e (0) = 48 K. The mask wavelengths are (a) 2 mm, (b) 
1.33 mm,(c) 1.0 mm, and (d) 0.5 mm. 



To perform spatially resolved spectroscopy, which pro- 
vides information on the ion velocity distribution, we spa- 
tially integrate the fluorescence signal given by Eq. [7]over 
the region of interest, 



dxdy F(v, x, y) 



region 



(11) 



(12) 



A particularly useful expression is found if the region 
extends over the entire x — y plane to collect signal from 
all ions illuminated by the laser sheet 
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where v x and u x are the ^-components of the veloci- 
ties, and we have evaluated integrals over v y and v 2 . 
If the hydrodynamic velocity arises purely from self- 
similar expansion (Eq. [5]), the spectrum reduces to a 
Voigt profifo 3 -^, 

loheff 



S p lasma{v) OC / dv 
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The rms width of the Gaussian component of this pro- 
file arising from Doppler broadening reflects both ther- 
mal ion motion and directed expansion and is given by 
s v exp — ^bTi + 7 2 c 2 , which increases with time because 
of plasma expansion (Eq. [SJ . 
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For a general hydrodynamic velocity, the spectrum 
may not explicitly take the form of a Voigt profile, but 
one can show that the width of the spectrum can be 
quantitatively related to the Lorentzian width, j e f f , and 
the mean square x-component of the velocity for all ions 
in the excitation volume, (i^) = s 2 /\ 2 + u 2 /X 2 . 

= 4aj 2 ff + (v 2 )/\ 2 , (15) 

where the overbar indicates a spatial average and the 
mean square x-component of the hydrodynamic velocity 
is 

ui= ( d 3 r^ I! ,(r) Ul (r) 2 . (16) 

The rms width of a Lorentzian is not well-defined, so the 
factor a depends upon the integration limits in Eq. 1151 
But a is constant and on the order of unity in our experi- 
ment, and the Lorentzian contribution to the linewidth is 
small compared to Doppler broadening, so any error in- 
troduced by this effect is small. When the hydrodynamic 
velocity deviates only slightly from the expansion veloc- 
ity, which is the case for small amplitude IAWs excited 
in this study, fitting the spectrum to a Voigt profile pro- 
vides good approximations of both the Lorentzian width 
and the the mean square x-component of the total hy- 
drodynamic velocity. 



C. Creating density modulations to excite ion acoustic 
waves 

In the plasma creation process, uniform intensities of 
both photoionizing beams, over the length scale of the 
MOT, result in plasmas with the same density profile as 
the MOT. However, we can spatially modulate the inten- 
sity of either ionizing laser and tailor the initial plasma 
density distribution. For these studies we have modu- 
lated the intensity of the second ionizing beam (412 nm) 
by placing a transmission mask along its path (Figure [2]). 
If we place the transmission mask before a single pass of 
the ionizing beam through the experimental region, we 
create high amplitude modulations in the plasma that 
access non-linear effects. For this study, however, we 
probe the linear regime by allowing the 412 nm ionizing 
beam to first pass through the atoms, then through a 
transmission mask and back onto the plasma. This cre- 
ates ~ 10% plasma density modulation with wavelength 
set by the period (Ao) of the mask. We translate the 
mask pattern to align a density minimum to the center of 
the plasma, and for this study, we used one-dimensional, 
50% duty cycle, binary patterns with wavelengths of 0.50, 
1.00, 1.33 or 2.00 mm (spatial frequencies of 2, 1, 0.75, 
and 0.5 cy/mm). 
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FIG. 4. ID slices through density profiles for T e (0) = 48 K 
showing density perturbations. The mask wavelengths are (a) 
2 mm, (b) 1.33 mm,(c) 1.0 mm, and (d) 0.5 mm. Dots are data 
and deviations from the fit Gaussian (solid line) represent the 
IAW density modulation, 8n (Eq. [T7J). 



III. RAW DATA SHOWING IAWS 

Figure [3] shows surface plots of the plasma density cre- 
ated with various periodic square wave masks and an ini- 
tial electron temperature of T e (0) = 48 K. We can clearly 
notice the modulations in density. Due to optical ac- 
cess limitations, the perturbations appear rotated with 
respect to the field of view of the camera (FOV), i.e. the 
ionizing beam and fluorescence beam intersect at a 74° 
angle in the imaging plane, so the imaging beam is 16° 
from normal to the IAW propagation axis. 

To study the evolution of the density modulations, we 
form ID density profiles from a central strip of width a by 
averaging the 2D data parallel to the modulations (Fig. 

By fitting such data to a Gaussian, we can separate 
the total density into background ncauss and wave 5n 
components 

rn(x, t) = n Gauss (x, t) + 5n(x, t). (17) 

Figure [5] shows the evolution of 8n{x,t) for a one trans- 
mission mask period with different initial electron tem- 
peratures. Note the oscillation of the wave in time and 
space. 



IV. EFFECT OF PLASMA EXPANSION ON THE IAW 

To motivate the formalism for analyzing IAWs, it is 
necessary to understand the effect of the plasma ex- 
pansion on the IAW amplitude. This can be found by 
treating the wave action as an adiabatic invariant 51 . We 
start by describing an ion-acoustic wave of wave vector 
k and amplitude 5n in an infinite homogeneous plasma, 
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FIG. 5. Evolution of 8n for mask wavelength of 2 mm. Time 
since ionization is indicated on the right, and Sn(x, t) has been 
scaled by instantaneous peak density naauss(0,t) and offset 
for clarity, for (a) T e (0) = 48 K, (b) T e (0) = 105 K. Solid red 
lines are fits to Eq. (|22|) . and the dotted black lines follow 
points of constant spatial phase of the standing wave. As 
the electron temperature increases, the frequency of the wave 
increases and for a fixed electron temperature, as the number 
of cycles per millimeter increases, the frequency of the wave 
also increases. 



for which the dielectric function is given by 
This implies that the wave energy density is 



(18) 



W k = — [ue{k,uj)\<—— = 19 
oui 8tt uo z 8n 

for wave field amplitude Ek, uj p i 3> w, and 2n/k <C Arj. 
The field amplitude can be related to the amplitude of 
the IAW through 
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The action density is given by Wj~/uj, which yields the 
total action in a volume V 
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For a slow expansion (\u/uj 2 \ <C 1) and if damping ef- 
fects such as Landau damping are negligible, the IAW 
will evolve adiabatically during the expansion, and the 
total IAW action will remain constant. For our mea- 
surements, |w/w 2 | <C 1 is a good approximation at early 
times, although for longer wavelengths and long times, 
\uj/ui 2 \ « 1. 

In our inhomogeneous UNPs, the length scale for 
density variations is large compared to the character- 
istic scale of the IAW (a 3> 2w/k), which suggests 
that Eq. [21] should apply locally^. This allows us to 
identify V(t) oc <r(t) 3 and oc l/a 3 (t). As shown 
previously^ and discussed below, k(t) oc l/a(t) and 



uj(t) » k(t)y/T e (t) I mi oc l/er(i) 2 , where we have used 
the time evolution of the electron temperature (Eq. [6]) . 
For constant Nk, this yields the final result for the evo- 
lution of the IAW amplitude \8n(x,t)\ oc l/a 3 (t) oc 
n-Gauss (0,f). So we expect the wave amplitude to scale 
with the the peak ion density, riGausss(0, t), and devia- 
tion from this behavior would indicate damping or gain 
processes. 



V. MEASURING IAW WAVELENGTH, FREQUENCY, 
AND DISPERSION 

To analyze data for IAW evolution such as to mea- 
sure the wavelength, frequency, dispersion, and damping 
rates, we fit the data to a standing wave model. As dis- 
cussed in the previous section, we scale the amplitude 
of the perturbations by the instantaneous peak density 
obtained from the Gaussian fit (Fig. [4]), and assume a 
damped standing wave with a Gaussian envelope, 



5n(x, t) 
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= A{t)e- x2 l 2 °°™^ cos [k(t)x 



cos [&(t)x] 



(22) 



This expression fits the modulation at a single time, t, 
and yields the instantaneous amplitude A(t), envelope 
size a e nv(t) and wavevector k(t). The hydrodynamic de- 
scription for an infinite homogeneous medium used to 
derive the IAW dispersion relation, Eq. [T] allows for pla- 
nar standing-wave solutions, but finite size, plasma ex- 
pansion, and density inhomogeneity introduce additional 
factors that are not small here and preclude an analytic 
solution, so the model behind Eq. [22] is only phenomeno- 
logical. The solid red lines in Figure [5] are fits to this 
evolution model in which A, a env , and k are allowed to 
vary independently for each curve. 

The dotted black lines in Figure [5] follow trajecto- 
ries of constant k(t)x. Notice that as the wave evolves, 
the wavelength/ wavevector appears to increase/decrease 
with time. If we hypothesize that the wavelength scales 
with plasma expansion (Eq. [3]) , then we have: 
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(23) 



Fits to k(t) measurements, with ko and T exp as fit param- 
eters, yield values for 2ir/kfit,o that match the transmis- 
sion mask wavelength. We also find very good agreement 
between the extracted fit value for r exp and the value ex- 
pected from self-similar expansion, Eq. [2j 

To show this behavior, we compare the evolution of the 
wavelength with the evolution of the size of the plasma 
by plotting, in Fig. [6] A(i)/Ao, ko/k(t) and a(t)/ao, 
where all quantities have been normalized to the values 
at t = 0. Initial wavelengths match the period of the 
mask used and all the data follow one universal curve, 
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FIG. 6. Evolution of the IAW wavelength and plasma size, 
normalized to initial values, for <jq = 1.45 mm and for (a) 
T e (0) = 25 K, (b) T e (0) = 48 K, (c) T e (0) = 70 K, and (d) 
T e (0) = 105 K. Size and wavelength follow a universal curve 
with no fitting parameters. For the solid line, r exp has been 
set to its theoretical value (Eq. [2j. The universal behavior 
for wavelength and plasma size indicate the wave is pinned to 
the expanding plasma. 

1/2 

(l + £ 2 / T eL P ) , with no fit parameters and r exp as cal- 
culated from <7o and T e (0) (Eq. [21). This universal be- 
havior for wavelength and plasma size indicate the wave 
is pinned to the expanding density distribution. 

To obtain a dispersion relation, we now need to ex- 
tract the frequency of the wave, ui(t). We can extract 
the frequency by analyzing the evolution of the pertur- 
bation amplitude. The amplitude variation in time can 
be modeled as 

A(t) = A e~ rt cos [<f)(t)] = A c- Ft cos (^J oj(t')dt'^j , 

where the integral accounts for accumulation in the 
phase. By fitting the evolution of the phase of the en- 
tire wave, we are essentially finding the frequency in the 
frame moving with each plasma element, and our mea- 
surement does not need to be corrected for Doppler shifts. 
For the frequency, we now assume the form of the dis- 
persion for an infinite, homogeneous medium from Eq. [I] 
the observed variation of wave- vector k(t) (Eq. I2U1) . and 
electron temperature evolution T e (t) predicted for a self 
similar expansion (Eq. [6]) to obtain 

w(t) = k{t)^k B T e {t)/m t = ( 1 + t l /T 2 ) • (25) 

Initial amplitude Aq, frequency uiq-, and damping rate T 
are allowed to vary in the fits, while r exp is taken as the 
value predicted for the electron temperature set by the 
laser and initial plasma size, (Eq. [5]) resulting in fits that 
match the data very well (Fig. [7]). As can be seen from 
Fig. [3 and Eq. [551 the frequency decreases with time. 
Similar behavior was observed in simulations of spherical 
IAWs±£. 



FIG. 7. Evolution of the amplitude of Sn(x, t)/riGauss(0, t) 
(Eq. [22]) and fits to obtain uj{t) for (a) T e (0) = 25 K, (b) 
T e (0) = 48 K, (c) T e (0) = 70 K, and (d) T e (0) = 105 K. Lines 
are fits to Eq. (|24[) in which r exp has been fixed to the theo- 
retical value. 
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FIG. 8. Dispersion relation of IAWs for different initial elec- 
tron temperatures. Lines are from the theoretical dispersion 
relation, Eq. [T] with no fit parameters. koXo < 0.42 for all 
conditions used in this study, with the average density being 
used to calculate Ad. 

We extract ujq and fcrj f° r a range of initial electron tem- 
peratures and mask periods, and calculate the dispersion 
of the excitations, as shown in Fig. [5] The excellent 
agreement with theory (Eq. [TJ) confirms that these exci- 
tations are IAWs. The planar standing- wave model cap- 
tures the dominant behavior of the wave, and to a high 
accuracy there is no deviation from the standard disper- 
sion relation in spite of the plasma's finite size, expansion, 
and inhomogeneous density. We suspect that this follows 
from the fact that the length scale for background den- 
sity variation and time scale for plasma expansion are 
reasonably long compared to the IAW wavelength and 
period respectively. Observations for longer times will 
be more sensitive to interactions between the wave and 
the boundary. 

In the first report of IAWs^, the frequency was found 
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FIG. 9. Evolution of the amplitude of Sn(x, t)/naausa(0, 0) 
for (a) T e (0) = 25K, (b) T e (0) = 48 K, (c) T e (0) = 70 K, and 
(d) Te(0) = 105 K, as was done for analysis in the first report 
of IAWs in UNPs£. 



by fitting the evolution of 8n(x, t)/ncauss{0, 0) (Fig. [9]) 
instead of Sn(x, t)/nG au ss(0,i) (Eq. [22]). Within our un- 
certainties, there is no difference in the resulting disper- 
sion relation. 



VI. IAW DAMPING 

The scaling of the amplitude has a significant effect on 
the interpretation of the damping of IAWs. The unsealed 
amplitude of the perturbations, 5n, decreases rapidly in 
time (Fig. [5]), but as discussed in Sec. IIV1 a significant 
decrease is expected due to plasma expansion, which is 
distinct from damping. This effect is accounted for by 
scaling the amplitude by the instantaneous peak density, 
5n(x,t)/nGauss{Q,t), which yields A(t) shown in Fig. [7] 
With proper scaling, it is clear that damping is small on 
the timescale of our measurements. 

The inverse of the damping rate, found by fitting Eq. 
[24l to the data, is shown in Fig. [TO] The typical damp- 
ing times are 3-10 times longer than the timescale of our 
measurements, and thus poorly determined. The obser- 
vation times were limited by decreasing density and sig- 
nal as the plasma expanded. Shorter wavelengths 4 ^ will 
increase the oscillation frequency and allow us to observe 
more periods and explore the damping with greater pre- 
cision. 

The damping of collective effects is interesting in plas- 
mas because it can be sensitive to collisional or kinetic 
effects^. Damping can also probe important many-body 
properties such as viscosity, which has become a topic 
of great interest in strongly-coupled systemsi^. In the 
case of IAWs, Landau damping is often the dominant 
damping mechanism. When T e ks Tj there is a large 
population of electrons traveling at just below the phase 
velocity of the wave that can efficiently extract energy 
from the wave. The Landau damping rate scales with 
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FIG. 10. Measured damping times (F _1 ) found from fitting 
the wave amplitude evolution to Eq. [24] Initial IAW wave- 
lengths are given in the legend. 



the oscillation frequency, uj (Eq. [T| , as given by^ 
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1L 



(l + fc 2 A^) 3/2 



D) 
1/2 



3/2 



exp 



-T e /Ti 



2(l + fc 2 Af 



DJ 



(26) 



To account for the time dependence of the damping rate 
over the evolution time of the experiment we average the 
calculated rate over the observation time, 



lavg = — I 1L (t) dt. 

to Jo 



(27) 



The inverse of this average rate gives theoretical damp- 
ing times of larger than 1 ms for all conditions studied 
here, which is much longer than the timescale observed 
for damping of IAWs. 

Intuitively, one expects other effects to lead to damping 
or dephasing of the wave in UCPs. For example, propaga- 
tion into low density regions at the plasma outer bound- 
ary should lead to decay. Also, when the wavelength de- 
creases and the dispersion relation deviates from constant 
phase velocity due to the nonvanishing Debye screening 
length (k\o rs 1), the dispersion relation should depend 
on plasma density and vary through the plasma. This 
will lead to a more complex wave evolution. 

It would be very interesting to access a regime in which 
Landau damping dominates. It is difficult to achieve sig- 
nificantly lower T e /Ti in UNPs because of intrinsic elec- 
tron heating effects at low temperature, such as three 
body recombination 34 . But at short wavelength, as k\o 
approaches unity, the damping increases sharply and 
should be observable. Much more work remains to be 
done to understand damping of IAWs in UNPs, nonethe- 
less, this work represents a valuable first step. 
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VII. ION VELOCITIES IN THE IAWS 

In the previous sections, we have studied IAWs through 
the variation of density in space and time. It is also 
possible to measure the velocity of the ions in the plasma, 
which should also display the effects of the oscillation. 

For total density given by Eq. [T7] with perturbation 
dn given by the standing wave model (Eq. I22j) , the ion 
flux J can be found from the one-dimensional continuity 
equation, dn/dt = —dJ/dx. For IAW wavelength much 
less than the characteristic size of the plasma (A <C a) 
and slow variation of a, ui, and A, the flux is 

J w ^-Ao72 Gauss (0,t)e~ x2/2cr ^e~ rt smkxsmu}t. (28) 
k 

The hydrodynamic ion velocity at x and t due to an IAW 
is given by uiaw — J/n. We observe that the envelope 
size is approximately equal to the characteristic plasma 
size ((Tenv ~ a). So exponential factors in this ratio can- 
cel, and 

uiaw ~ -7-A e~ rt sinfcx sincji, (29) 
k 

where we have assumed a small perturbation. In addi- 
tion to the IAW velocity, ions possess the hydrodynamic 
velocity due to plasma expansion (Eq. [5]) and random 
velocities characterized by the ion temperature Ti. 

Section III Bl describes how the Doppler-sensitive laser- 
induced fluorescence spectrum can be used to obtain the 
mean square of the velocity along the laser propagation 
direction for all ions in the excitation volume, (v%} (Eq. 
IT51) , where the overbar indicates a spatial average. For a 
small-amplitude IAW in an expanding plasma, («§) arises 
from thermal, expansion, and IAW contributions, 

(v%) = k B Ti/mi + -i 2 a 2 + u 2 IAW + 2jxu IAW . (30) 

The average of the cross term, 2^/xuiaw, tends to zero 
for small IAW wavelength, and is already small for our 
conditions. We will neglect it in subsequent analysis. 
For clarity, we have also neglected effects of the small 
angular misalignment of the IAW and laser propagation 

directions. 

Figure [TT] shows measurements of (v%.) extracted from 
the spectra for T e (0) = 70 K expressed in terms of an 
average kinetic energy of the ions, mi(u|)/2. The small 
offset for the kinetic energy at early times arises from 
disorder-induced heating of the ion o 32 i 36 . The kinetic en- 
ergy then shows a gradual increase on the timescale of 
r exp as electron energy is transferred into expansion ve- 
locity of ions. Finally, it approaches a terminal velocity 
~ fcsT e o/2mi when all the energy has been transferred. 
Note that the overall evolution of kinetic energy is af- 
fected very little by the presence of IAWs, which further 
confirms that the plasma expansion is not affected signif- 
icantly by the small IAWs. The only deviations are small 
oscillations at early times. 
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10 20 30 40 

Time [lis] 

FIG. 11. Evolution of the average ion kinetic energy for 
T e (0) = 70 K, and various IAW wavelengths or no IAW. Data 
are fit to — ffcsli + \m,i^ 2 a 2 . This emphasizes that 

the data follows a universal curve indicating a negligible effect 
of the modulations on the overall expansion. 

Figure [T2] sh ows the difference between the data and 
the fits to \mi{v1) — \ksTi + \mi^ 2 a 2 , which empha- 
sizes the effect of the IAW. This difference is now fit to 



AE KE = -m,iU 2 IAW = -m, [T A ° e J sin2 ^ ( 31 ) 

plus an arbitrary offset, where the evolution of u/k is 
replaced by the expected IAW speed \/ KbTJ^) / 'mi, the 
evolution of u is set by Eq.[25l and the only fit parameters 
are the fractional amplitude Aq and decay rate T, Notice 
the clear oscillations in the kinetic energy at twice u> that 
are out of phase with the oscillations in amplitude of the 
density perturbation shown in Fig. [7] 




10 20 30 

Time [lis] 

FIG. 12. Oscillations in the ion kinetic energy for T e (0) = 
70 K, extracted from the data shown in Fig. 1111 obtained 
upon subtraction of the calculated average ion kinetic energy 
from the fits shown in Fig. 1311 The solid lines represent a fit 
to Eq. [31] and for clarity, the curves are offset. 

Figure [T^] can be interpreted as a measurement of the 
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amount of energy in the waves. This energy oscillates 
back and forth between potential energy of ions in the 
electric field when the density perturbation is largest and 
kinetic energy when the kinetic energy is largest. Mea- 
suring velocity perturbations may prove to be a more 
valuable probe of IAWs than density perturbations when 
the wavelength is very small and difficult to resolve opti- 
cally . 



VIII. CONCLUSIONS 

We have presented expanded studies of ion acoustic 
waves in ultracold neutral plasmas- , including a descrip- 
tion of the effects of plasma expansion on the ampli- 
tude of the wave. We observe the wave through density 
and velocity perturbations, and we find excellent agree- 
ment with the well known IAW dispersion relation. The 
IAWs damp significantly faster than expected for Lan- 
dau damping, and this effect remains to be explained. In 
the future, we plan to extend these measurements to the 
short wavelength 49 , non-linear dispersion regime of ion 
plasma oscillations, where strong coupling is predicted 
to have an effect and Landau damping increases. With 
these tools we can also study different initial plasma ge- 
ometries to access other phenomena such as streaming 
plasmas, nonlinear waves, and possibly shock physics. 
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